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We calculate and explore the moduli potential for M theory compactified on G 2 -manifolds 
in which the superpotential is dominated by a single membrane instanton term plus one 
from an asymptotically free hidden sector gauge interaction. We show that all moduli can 
be stabilized and that hidden sector gauge symmetries can be Higgsed at a high scale. We 
then compute the spectrum of superpartner masses at the GUT scale and evolve it to the 
electroweak scale. We find a spectrum which is very similar to the G2-MSSM with light 
gauginos - accessible at the LHC - and a neutral wino dark matter candidate. 
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I. INTRODUCTION 

String/M theory vacua with four large spacetime dimensions are effectively described by a low en- 
ergy theory which typically has non-Abelian gauge symmetry, chiral fermions and hierarchical Yukawa 
couplings. These facts were essentially evident in 1984 p]. Since these are some of the key ingredients 
of the Standard Model of particle physics, this provides some essential motivation for pursuing questions 
in so-called " string phenomenology" . 

At that time, the vacua mainly under consideration were exactly supersymmetric - a condition moti- 
vated by the relevance of supersymmetry to the stability of the mass scale of the Standard Model under 
quantum corrections. The problem of how supersymmetry is broken is closely related to the moduli 
problem(s). Supersymmetric string/M theory vacua have moduli fields - namely, the extra dimensions 
have zero modes in the supersymmetric limit. These zero modes are continuous parameters of the metric 
on the extra dimensions (e.g. Calabi-Yau or G^-holonomy metrics). They parametrize the size and shape 
of the extra dimensions and appear in the low energy theory as fields with Planck suppressed interactions 
with matter. Having massless scalars with Newtonian couplings is problematic phenomenologically for a 
number of reasons such as "fifth force experiments" and cosmological constraints. 

The standard lore was that, once supersymmetry breaking is "understood" , these moduli would receive 
masses from supersymmetry breaking effects. Since all couplings and masses are functions of the moduli 
fields in string/M theory any supersymmetry breaking effects in the effective Lagrangian will depend 
upon them. Therefore, whatever mechanism is responsible for supersymmetry breaking will necessarily 
also generate a contribution to the moduli potential. If supersymmetry is the solution to the hierarchy 
problem, then such contributions to the moduli potential must have a size of order M^ sy (where the 
value of M susy depends upon how supersymmetry breaking is "mediated" to the observable sector). In 
any realistic case, M susy is orders of magnitude smaller than, say, the Planck scale. One of the simplest 
ideas for generating a small scale of supersymmetry breaking, which goes back to Witten [2] and others 
is that this small scale is the strong dynamics scale of a new hidden sector interaction, in close analogy 
with the dynamically generated scale of QCD. 

Putting all of these ideas together, strong dynamics in a new hidden sector interaction should i) 
generate the hierarchy between m p i and M\y, ii) break supersymmetry and iii) generate a potential 
which stabilizes all the moduli fields. Attempts at realizing these ideas in string theory have met with 
partial success over the years. Nilles [3] made the earliest investigations in the supergravity context, [HE] 
considered the heterotic string. In all of the investigations in string theory, the problem could only be 
considered to be partly solved in the sense that it was not understood how to stabilize ALL of the moduli 
fields whilst simultaneously generating the mass hierarchy. 

The problem was finally given a fairly satisfactory solution [H] in the context of M theory vacua in 
which the seven extra dimensions form a manifold of GVholonomy [TJ [8] . There is a simple explanation 
for this success: in M theory, in the absence of fluxes, all the moduli fields reside in superfields which 
undergo shift symmetries (because they all contain axions). This implies that only " non-perturbative" 
contributions to the potential are allowed - such as strong dynamics of the type we are interested in for 
the hierarchy problem. Thus, the potential energy will always be exponentially small in M theory. This 
is not the case in the heterotic string and other string theories. For instance, there are no symmetries 
preventing perturbative contributions to the potential for complex structure moduli in the heterotic or 
Type IIB string theories. 

In |6j the main idea was to consider hidden sectors which consist of two or more asymptotically free 
gauge sector^] The purpose of this paper is to consider the effect and role of other types of terms in 
the superpotential on the moduli potential. In M theory all terms in the potential are non-perturbative. 



1 One hidden sector gauge group is actually enough to stabilize all the moduli, however the vacuum is not in a region where 
the supergravity approximation is valid. 
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For instance an ordinary cubic Yukawa coupling between three charged chiral superfields is generated by 
a membrane instanton which wraps a 3-dimensional submanifold of the extra dimensions which passes 
through the points where the fields are localized [18J. As such it will be exponentially suppressed by 
the volume of this cycle. We will consider a hidden sector U(l) gauge group with charged matter with 
a cubic Yukawa coupling as well as a single asymptotically free hidden sector gauge group (which turns 
out to be necessary for obtaining vacua treatable within the supergravity approximation). 

We will follow the ideas of [6] to understand the minima of the supergravity potential. We find that 
all moduli can indeed be stabilized by these contributions. We then go on to calculate the observable 
sector supersymmetric partner masses and couplings at the electroweak scale. We will see that this 
extends the framework of [6] in the sense that it permits for a wider range of gravitino masses (i.e. 
supersymmetry breaking scales), however the low energy phenomenological features for observationally 
consistent GVmanifolds remain essentially the same as those in the [6] . The discussion section at the end 
of the paper briefly summarizes the main phenomenological features and suggests some further directions 
to investigate. 

II. FOUR DIMENSIONAL EFFECTIVE THEORY AND MODULI STABILIZATION 

The complexified moduli space A4b 3 of M theory compactified on a real manifold manifold Xj of G2 
holonomy is parametrized by holomorphic coordinates 

Zi = U + isi , (1) 

where i runs from 1 to 63, the 3 rd Betti number of Xj. Moduli and axion fields, respectively denoted by 
Si and ti, arise from the 11 dimensional metric and 3- form field as Kaluza-Klein zero modes of Xj. These 
fields pair up to form the bosonic components of chiral superfields in four dimensions^] All the parameters 
of the low scale physics are functions of the moduli vacuum expectation values, so stabilizing the moduli 
is a fundamental step towards phenomenological implications. As an side, note that the axions, being 
periodic fields, are not moduli. 

A. Moduli Stabilization and Symmetry Breaking 

The tree level Kahler potential of the effective four dimensional supergravity theory is [9] 

K = -31n(4^ 1/3 y x ) , (2) 

where Vx is the volume of X7 in units of the eleven dimensional Planck length. Vx is a homogeneous 
function of moduli of degree 7/3. It has been demonstrated in |10j that the observable sector matter 
spectrum (e.g. superpartner masses) depends very weakly on the detailed form of the moduli Kahler 
potential - a fact due to the homogeneity of Vx- 

We will consider Gi manifolds which give rise to i) a hidden sector gauge theory which confines at 
low energies ii) a hidden sector U(l) gauge group under which there are three charged matter superfields 
T a with a cubic superpotential. The confining gauge theory arises from a co-dimension four orbifold 
singularity along a 3-dimensional sub- manifold of X-j. Due to the axionic shift symmetry and holomorphy, 
the coupling of the cubic term in W is generated by membrane instantons and is therefore exponentially 
suppressed by moduli vacuum expectation values. 



2 The zero modes of the eleven dimensional gravitino field along X7 gives rise to the four dimensional gravitino and 
modulino-axino fields required by supersymmetry. 
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The total superpotential is thus given by 

W = A e 27Tif/p + A TiT 2 T 3 e 27Tif ~ . (3) 

In the above A (due to PQ shift symmetry) and A are 0(1) constants. The constant P is the beta function 
coefficient of the confining gauge theory. These two terms are regarded as the leading contributions to 
W. The superpotential can, in principle contain many additional non-perturbative contributions if X-j 
contains additional supersymmetric 3-cycles. 

The gauge kinetic function of the gauge theory in the hidden sector is a linear combination of the 
moduli and given by 

b x 

/ = E^ = ^+^5s- (4) 

8=1 

The imaginary part is the inverse of the hidden sector gauge coupling a^g which plays a central role in 
this work. Since, in the UV hidden sector is weakly coupled, we can regard the coupling as an expansion 
parameter in the eventual vacuum we will study. This is consistent as long as the vacuum has moduli 
vacuum expectation values such that ckhs C 1. / is also a linear combination of moduli 

"x 



~f = Y J N i z l , (5) 



i=l 

the imaginary part of which is the volume of the 3-cycle wrapped by the membrane instanton. iVj and iVj 
are positive integers which are determined by the homology classes of the 3-cycles. These are microscopic 
parameters determined by the details of the compactification manifolds. 
The Kahler potential for the matter fields in the U(l) sector is given by 

3 rri rrs 

rUi = y ±a±a (g) 

^ V X 

a=l 

This moduli dependent form has been proposed in different ways for calculating the moduli dependence 
of the matter kinetic terms |10j. 

The scalar potential in 4-dimenasional N = 1 SUGRA is determined by 

V = e G (j2G I G I -3) . (7) 

The summation is over all the moduli and matter fields (/ = i, a). The Kahler function G is defined as 

G = K + ln\W\ 2 . (8) 
The F-terms are given by Kahler covariant derivatives 

Gl = ^r = ^r + dlK ' (9) 

and G 1 = J2j kIJg j where rIJ is the 

inverse Kahler metric. 

Whereas there are two terms in the superpotential, as a cross term, only one linear combination of 
axions ti and phases of the complex scalar fields r a appears in the scalar potential. Minimizing the scalar 
potential fixes this combination to be (for A > 0) 

cos 



WiNiti - Ti - T 2 - 73) = -1 , (10) 
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where we have defined tOj as 

Wi = — . (11) 

In order to stabilize the rest of the axions, we need to consider some 63 + 2 more subleading terms 
in the superpotential generated by membrane instantons (obviously with larger actions than the term 
considered above). The extra cross terms generate a potential stabilizing the remaining axions and is 
studied in detail in [TT] . 

In order to simplify the study we introduce the following notation 

x=jT 01 T 02 T 03 e^ N ^ . (12) 

The superpotential can be rewritten more compactly as 

W = A(1- x) e-KN-s/P^t/P (13) 

The F-terms can be derived explicitly as follows 

2vr 2vr 

Vi = — — x — . (16) 

In the above and in the following, there is no summation over repeated indices i and a unless explicitly 
written. We have also introduced some variables a\ defined by 

_ ainVx n7 , 

such that en = 7/3. We then can write 

dK _ 1 dK _ _ 3 d\nV x = _ 3 en ^ 

dzi 2% dsi 2i dsi 2i Si 

Using the Kahler potential Q, one can obtain the components of the Kahler metric. Using the 
orthonormality condition 



where 



Y,K I3 K JK = 5 I K , (19) 



one can read the non-zero components of the inverse Kahler metric. After some work, one obtains the 
following expression for the raised index i^-terms 



G i = _ iSi ( 4 1 y J^ a + 2— + 2] 

KZai^l-X 3 JV ' l-X J, , E,T.Tq 

n 1 1 qi/„ 



(20) 

a 



G T a = _ T ( 4 1 x X X I 2 1 - V m c N 7 1 TaTa I 7 1 X 
a \3 l-xT a T a 3i 1 HgTaTg i-x 1 1 9i 1 S„7Vr B U x 9, , S.t.t, 1 - g 



(21) 
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In the above, the matrix Aj,- is defined as 



1 



1 /l 



(22) 



where the matrix Kij is the Hessian of K which is related to Kahler metric by a factor of 4. One can 
easily check that the matrix Ay satisfies the following properties 



"x 



b 3 x 



£ ^ = 1 

i,j =1 



3=1 



a; 



(23) 



Finally, putting all the pieces together one obtains the following expression for the scalar potential 

\ r r Z_* ... I 



pEo T a T a IVx , „ 

F(,, t ,r Q ) = ^ 2 (l -x) 2 e-**' p X<*'< 



3 1 , E a ?qr a aj 1 — £ 1 

3Vx 1 J 

+4— -= V NiSi- 

E 



(1-x) 2 



Vy 7 1_ 



1 — x 



14 



1 ^-v T a T n 



9 ! , EJ a r» ^ V x 

3Vx a 



4yT a r a 7 

3^ F X 9i | E a ^qTq 
a 1 3Vx 



1 / T a T a 



+ 4 



(24) 



In order to study the vacua of the potential, we make the following ansatz for the moduli vacuum 
expectation values at the minimum 



Cli 1 — x 
Si = — s . 



(25) 



The unknowns aj and s will be determined later. Then, the inverse coupling of the hidden sector gauge 
theory is given by 



*> 3 x 



i=l 



a, 



(26) 



In order to study the minima of the potential, we use the fact that the hidden sector gauge theory in the 
UV is weakly coupled 

aus « 1 • (27) 

We keep working in this regime and regard the gauge coupling «hs as our expansion parameter and 
expand our solutions in its powers. We understand from (26) that the weak coupling limit is attained 

(28) 



when one or more of the j/j, 1 < i < b x , are small 



Vi ~ «HS 0, i G x ■ 
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For the moment, we also assume that s is non-zero and finite in this limit, a crucial assumption to be 
verified later. For simplicity, we use the following parametrization 



= -Ci«HS = -C«HS + (C - Ci)«HS, i^X, 



where £ and Q are 0(1) constants. From the definition of yi (16) and by the following expansion 



Ni 



we find that 



and 



= n + njQHs + 0(«hs)) iG X 



2vr 2vr 
y = -C«hs = — - x 



P 



n 



C-Ci 



27r nj 
n 



(29) 



(30) 



(31) 



(32) 



The constant £ is a measure of the relative difference between the two homology classes which generate 
the two terms in the superpotential. One can show that it is less than 7/2 (for the case all of the yi go 
to zero). Furthermore, in this parametrization we can also write Wi as follows: 



2tt 2tt P 

Wi=w + u^ohs + C(aHs) = S h (C ~ Ci)~ «HS + 0(«hs)) 

P n n 



Finally, from ( 16 ) we can show that 



x=£(l + C=|^)+0(4.)- 



The gauge coupling in the hidden sector can be re-expressed as 



a 



-l 

HS 



1-X 



c»£? + (i-*) 5 £^i^cs£ 



2/ ^ ^ — ' Ci * y% y 

Thus, the ansatz for the moduli in the two sets can be rewritten as follows 



1 Ot/Ci -1 
a HS ' 



4 ^E^/O 

1 CLi/Vi _i 

Moreover, from ( 12 ) one can see that 

- ^WiNiSi = In (~jr T m T Q2 T( 

i ^ 

which to the leading order in ohs can be written as 

1 



(33) 
(34) 
(35) 

(36) 
(37) 



(38) 



-i 



P A 



(39) 
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We take this as a constraint. For positivity of the coupling, either of the two following conditions have 
to be met 



n n A 

i: — < 1 and — < -^ToiTo 2 To3 , 

n n A 

it : — > 1 and — > -jT 01 T 02 T 03 . 



(40) 
(41) 



In order to determine the moduli vev's, in the first step one must determine the ai. Putting back 
the ansatz (25) into the definition of aj (17) one gets a system of 63 transcendental equations which, in 
principle, completely fix the 



(l; 



d(-K/3) 



din Si 
8(-K/3) 



. 1 h/U n ,-i 



d In Sj 



_ 1 a i/Vj — 1 



,i i X ■ 



(42) 
(43) 



For the moduli to be positive, for positive iVj, the solutions for dj must all be positive. One can numerically 
check that for quite generic sets of parameters the above system of equations yield positive solutions for 

Next we determine s in the limit of interest, i.e. m —> 0, and verify the initial assumption. We 
minimize the scalar potential ([7]) with respect to the moduli Si 



dV_ 

dsi 



. 



(44) 



We multiply each of these ^3 equations by Sj and sum all of them together. This gives an equation for 
s which can be used to determine qhs an d this, in turn can determine the individual moduli vacuum 
expectation values from our ansatz (25) 



„ w 



C-xs z 



X 



3 y 



+ 4C 



„ ± a ± a 



W X 
S 

yl-x 
7 1 



1 I SZq TgTg 

1 + 3V X 



3 I I En TgTg J 1 

1 + 3V X 



14 



T a T a 

9 , , HgTgTg ^ V X 

1 + 3V X a 



E 



(45) 



Solving the above equation to the first subleading order in y results in 



2 1 -x 



+ 



l-x 



5 3 
14 + 14T 



E 



Vx 
T a T a 



1 + 



TgT a 

3V X 



y_ 



(46) 



We can see that this solution is non-zero and finite when y — > and so is self-consistent. This is the 
solution which describes the minimum of the potential. The second root which goes to zero in the limit 
y — > corresponds to a maximum of the potential. 
In the minimum, one reads the F-terms as follows 



Gi = iNiyi 



1 



1 + 



3V X 



1-n/P 

_T a f n/P / a HS -P \ Vx \ 



(47) 
(48) 
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and also 



G l 



-IS) 



10 2 n/P V x 



n/P ^aus-P 



2vr 



G 



T a T a 



3 1 + J^TaTal-n/P 21 71 - n/P ^ T a T a J I -n/P 

4 n/P V x _ 7 1 3 - n/P _ 7 1 
3 + l- n/fTJ,"9i | Y. a TaTa l-nlP~ vV7 j: a T a T a V x 

1 ^ 3V X 1 ^ 3V X 

n/P V x 5 1 7 5 1 n/P ^ F x 

' a \l- n/PT a T a ~9 1+ l _ n/P + 9 + 3 1 - n/P \ j\T a 

3V X 

7 n/P 1 ^ n/P ^a HS ■ P 



(49) 
(50) 



91-n/P-i , J2gTaT a ) l - n /p s 2vr 
1 + 3V X 



The expression for s in the minimum is 



3 1 /5 1 /5 3 n/P V Jx W Tj«T fl \\ 

21 -n/P V 21 -"/^ U4 141-n/P^r a T a /V / / 1 



n/P ^a H s--P 



n/P^ 2vr 



(51) 



The scalar potential at the minimum of the moduli to subleading order in ckhs is read as follows 

= ^l e J2 a TaTJV x -7 -4,/ aHS -P x 
I mm 647T 

1 



X 



4(1-^(1-^^)-^ 



3Vx 
T a T a 



v 1 H — « 



T a T a 



(52) 



- n/P) 2 (l - 2^-C^^) T 
3 V 1 ' V 1 -n/P s 2tt / ^ 

+ £<Wi?(i-A^) 



n/P A Q H S • \ - PgPa 

Vx 
1 



E 



n/P^ 2tt / , , T, a TaT a ^ \ V 



TaTa 



B. Supersymmetric Vacua 

We first consider the existence of supersymmetric vacua. The condition for a supersymmetric vacuum 
are vanishing P-terms 



Gi = , 
G Ta =0, 



which implies 



3 a,; 1 



1 1 ai/yi Pln (n i T oi7o2T 03 ) 



TaTa 
V X 



2 iVj Vl 2vr JVj a^/y* P/n - 1 

x n/P ( ^ i ^ans ■ P N 



1 - x 1 - n/P 



(53) 
(54) 



(55) 
(56) 
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From above, for supersymmetric solutions to exist, it is necessary to satisfy 

(n/P) < 1 . (57) 
For the moduli to be positive, the constants have to satisfy the conditions of the two different branches 



(40) or (41). Thus, the SUSY condition (57) is consistent with branch i. 



C. Metastable de Sitter Vacua 

Here we turn to non-supersymmetric minima. In order to find the vacua and the vev's of the T a , we 



minimize (52) 



dV\mm{T a T a /Vx, a-Rs) _ , . 

d{T a T a /V x ) 

Thus, we get a system of three equations 

(l n/P) 2 (l 2 aus ' P \( 4 7 1 T a T a f 1 T a T a 

"~ n/ ' \ 1 ~ 2 l-n/P 2tt J { ~ 3 + 27^ ~E a r a r a y x \] E a ^ Ty~ 

i+ w x i+ 3V X 



1 - n/P 2vr J 2777 E a r a T a I J2 a T a T a V x J ( 59 ) 

1 3V X 1 ^ 3V X 

, /p)2 ( 1 2 (xhs-P \Vx(Vx 7 1_ T a T a \ 7 1 = 

W j V 1 - 2vr T a T a \T a T a 9 j , E a T a r a y x ^ 3 / E„t„t/ 2 " 

^ 3Vx ^ 3Vx 

Furthermore, it should at least be possible to tune the cosmological constant to zero for phenomenological 
reasons 

V\min(T a Ta/Vx,aus) =0 , (60) 

This gives the condition 



(WP) y- 2 i-n, P 2. ){ 4 -- 3 2:^r+- 9 ^tX^) 

a J. T 3 y x a 

+ n P(l-n P)[l + — 8 = > -^1 (61) 

1 1 ' ' V 1 - n/P 2vr J V 9 I 1 J2 a T a T a ^ V x J { ' 

1 "•" 3Vx a 

+ <„/iWl + — V ^)(VJ^ + 7 1 \ , 1 



1 - n/P 2tt y V ^ T a T a 3 1 , E a ^ a r a J J2 a T a T a 

7 V a 1 + 3V X 7 1 + 3V X 

One can actually solve the above equations for T a T a /Vx and ans to find the condition for a solution 
in a de Sitter vacuum as a function of n/P. We have solved them both analytically and numerically. 
Meta-stable de Sitter vacua only exist for the case n/P < 1, not for n/P > 1. The solutions are given as 
follows 

^ = ^^(1 + (9/14 -n/P) /(n/P)) , (62) 
a H s-P 11 -n/P/ n/P /14 h, , , N 
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The functions f(n/P) and g{n/P) have relatively complicated, uninteresting forms, they are plotted in 
figure [l] We also plot TT/Vx and cyhs versus n/P in figure [2J To further simplify the future formulae 
we use approximations to the vacuum expectation values and cvhs as a function of n/P, which are also 
given in the figure. 




0.0 0.2 0.4 0.6 O.S 0.0 0.2 0.4 0.6 0.8 

n/P n/P 



FIG. 1. Functions f{n/P) (left) and g{n/P) (right) versus n/P. 




0.0 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

n/P n/P 

FIG. 2. T a T a /Vx vacuum expectation values (left) and «hs (right) versus n/P. The dashed lines approximate 
those functions which are good enough in the range 0.2 <n/P< 0.8. 



Independent Numerical Calculation We have numerically checked that our formulae give the correct 
results as a cross-check that our methodology is robust. We give here an outline of a simple set of 
examples so that the reader has an idea of the numbers that we obtain. 

We choose a compact G2 manifold with 2 moduli and volume 

Vx = j^c k f[sf , (64) 
fc=l i=l 

with the following microscopic parameters 

ci = 1, c 2 = 1/3, c 3 = 1, 

a\ = 7/6, a\ = 7/6, a\ = 1, a\ = 4/3, a\ = 1/3, a\ = 2, 
N 1 = N 2 = 1. 

We scanned over all parameter space. The vev's of the moduli and scalar fields for different values of 
n/P that we obtained are of order 

(si)« 10 - 50, (65) 
(TT/Vx) ~ 1.1 — 0.3 . (66) 
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There was good agreement with our formulae. The fine-tuning of cosmological constant is done through 
X/A which take values between 0.1 and 100. 



III. PHENOMENOLOGY 

In this section we will study some of the phenomenology of the G2 vacua with emphasis on the soft 
SUSY breaking parameters. The values of the soft SUSY breaking parameters at the unification scale is 
determined in the standard way; the moduli and the F-terms are replaced by their vacuum expectation 
values in the supergravity Lagrangian. We then take the limit m p i —> 00 while keeping rrc.3/2 fixed to obtain 
the soft supersymmetry breaking Lagrangian. Having determined soft parameters at the unification scale 
we will RG evolve them to the electroweak scale. 



A. Stabilized Volumes 

There are basically two volumes which control the low scale phenomenology, the volume of the com- 
pactification manifold Vx and the volume of associative 3 cycle V l3 . They are homogeneous functions 
of the moduli of degrees 7/3 and 1 respectively. Since we are interested in the weak coupling limit and 
V la = e%s> we ex P ec t that 



Kr-II- 



c a 



-7/3 

HS 



and 



V, 



73 



O 



HS 



(67) 



In the above c is a constant which plays a useful role [12]. The first relation implies that the linear 
combination of moduli which determine the volume of the 3-cycle 73 dominate the volume function. We 
have numerically checked that, for generic compactification manifolds, it is a good approximation for 
Vx- Once the moduli get vacuum expectation values the volumes also become fixed. Phenomenological 
parameters of the low energy physics are functions of these volumes. They are plotted in figure [3} 



40 - 




20 - 



1400 




FIG. 3. Volume of associative 3-cycle V l3 (left) and volume of the compactification Vx (right) versus n/P. 



The dimensionless parameter c in (67) is 



V x 



V 



7/3 



(68) 



73 
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and it is useful to have an estimate of its range of values. We combine the two relations in (67) in two 
extreme cases in order to get bounds on c. In one limit we take all Oj's equal to 7/3N and in the other 
limit we take one of a^'s, say oq, much bigger that the rest. It gives the following bounds 

0(1(T 2 ) ~ (NN)-^ < c < (|^) 7/3 ~ 0(1) , (69) 

where TV" is the number of moduli and N is the geometric mean of A^'s. Note that these bound are 
coming from microscopic properties of the compactification manifold X. 

Similar bounds comes from phenomenological requirements. Loosely speaking, validity of the super- 
gravity approximation requires 

Vx > 1. (70) 

On the other hand, in order to make M theory (KK threshold) corrections to the parameters of the 
four dimensional physics negligible at and below m-GUT) we demand 11 dimensional Planck mass mu to 
satisfy 

mu > m GUT ~ O(10 16 GeV). (71) 

From dimensional reduction, one finds the following relation between the 4 and 11 dimensional Planck 
mass (m pl = 2.43 x 10 18 GeV) 

™ P i_(v x \' (72) 



Thus, one finds the following upper bound on Vx 

V X <J^-Y . (73) 
Wgut/ 

All things together, we need Vx to satisfy the following constraint 

0(1) < Vx < O(10 4 ) . (74) 

Again, the lower bound comes from requiring validity of supergravity regime and the upper bound comes 
from requirement of having a high compactification scale. It is basically 

m pi Z m n > m G u T . (75) 

To make these phenomenological constrains consistent, one demands 

0(5) < V 13 « «gut £ 0(100) • (76) 

As an example, for the case of the 3-cycle supporting the Standard Model gauge group, one finds the actual 
value at which the volume Vx must be stabilized can be determined from the GUT scale («qut = ^73 ) 

n= vx__ L^fJ^LY M 7 « v in -^^W 2 - 2xlQl6GeV ^ 2 

,7/3 



^"gutM^H ~7.51xl0- 2 9^ <0(1)- (77) 

y^i* 4vr GUT U G ut/ V 0.04 J V m GU T / 
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B. Gravitino Mass 

The gravitino mass plays an important role in gravity mediated models of SUSY breaking and sets 
the typical mass scale in the supergravity Lagrangian. The bare gravitino mass is defined as 



K/(2m 2 ,) I t t t" / . 

m 3/2 = m pi ' e p J \W/m 



(78) 



Having tuned the tree-level vacuum energy, we can determine gravitino mass using 



e J2 a TaT a /2V x 

™ 3 /2 = m pl —^^r- A(l - x) e~ 2 ^ p 



0.035^c~ 3/2 (1 - n/P) (l + 0.46 "'* Pe' 

V 1 — n/P 



n / P p -4A(n/P)\ -ll/2(n/P)-7.66e 4 - 4 ( n / p )/P 



(79) 



m p i 



As can be seen, the gravitino mass is determined by four parameters, n/P, P, A and c. The following 
plots in figure shows the behavior of m^j2 versus n/P. 




P=10, AiT 3 ' z =100 




0.82 



FIG. 4. Gravitino mass versus n/P. The lower dashed line is 10 TeV for reference to the cosmological BBN bound 
(30 TeV) and the upper dashed line is from non-thermal moduli dominated cosmology (500 TeV). 



C. Embedding the MSSM 

The MSSM, being a gauge theory with charged chiral matter, resides on a 3-manifold (with orbifold 
and conical singularities) different from the one supporting the hidden sector. Since two 3-manifolds in 
a 7-dimensional manifold generically do not intersect, the supersymmetry breaking generated by strong 
gauge dynamics in the hidden sector is generically mediated to visible sector via Planck suppressed 
interactions i.e. it is gravity mediated. In the analysis in the following sections, it is assumed that there 
is a GUT gauge group in the visible sector which is broken to the SM gauge group by background gauge 
fields (Wilson lines), following Witten |13] . 

The full low energy supergravity theory of the visible and hidden sectors at the compactification scale 
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(unification scale) is defined by the following Kahler, superpotential and gauge kinetic functions 

K = K + K Ul + K vs = -3 In (Att^Vx) + + K^{ Si )^¥ + (Z( Si )H u H d + h.c.) , (80) 

W = W us + W Ul + W ys = A e 27Tif /p + A TiT 2 T 3 e 2wif ~ + Y^ s ^ a ^^, (81) 

b 3 b 3 b 3 

f = J2N^ Zl , f = ±N lZl , f vs = £ N ? s Zl , (82) 

i=l i=l i=l 

where $'s are chiral matter fields, K^p is the Kahler metric and Y a ^ are un-normalized Yukawa couplings 
in the visible sector. Evidently, there is no bilinear (mass) \i term in the superpotential. It can be 
forbidden by a discrete symmetry which ensures doublet-triplet splitting [13] . This term can be generated 
by non-renormalizable terms in the Kahler potential (via Giudice-Masiero mechanism [16]) or in the 
superpotential. 

The Kahler metric KYq for the visible sector matter is assumed to be approximately flavor-universal 
and given by 

This can be ensured by U(l) symmetries associated to the vanishing 2-cycles at the location of the matter 
multiplets [17] . Higher order flavor- non-diagonal contributions to the Kahler potential can come from 
gravitational couplings between the scalars T a and the matter fields $ Q in the visible sector, but these 
are expected to be small due to the U{1) symmetry. Hence, the Kahler metric of the visible sector is 
diagonal and the Kahler potential can be written as 

K YS = K^Sap® 61 ^. (84) 

The gauge kinetic function in the visible sector is an integer linear combination of the moduli fields. 
The integers N ys is fixed by homology class of the 3-cycle supporting the observable sector gauge and 
matter fields. 

The Yukawa couplings arise from membrane instantons which connect singularities where chiral su- 
perfields are supported and given by 

Y aPl = c aM e 2m ^ 1 ^ , (85) 

where c ~ 0{1) and l % G 7L. Thus, due to exponential dependence on the moduli it is natural to obtain a 
hierarchical structure of Yukawa couplings [18l [19] . 

For simplicity, here we assume an observable sector which is precisely the MSSM, though GVmanifolds 
with additional matter can be treated in exactly the same way. Having all the moduli fixed, the F-terms 
can be explicitly computed in terms of the constants which determine the moduli potential e.g. n/P. 
Given a particular G2-manifold and a particular set of microscopic constants one can obtain a particular 
point in the MSSM parameter space. This defines a particular "G2-MSSM" . 

D. Physics at the Unification Scale 

1. Gaugino Masses 



In general three different sources contribute to the gaugino masses mi/2, the tree-level and anomaly 
mediation contribution as well as threshold corrections. They can be either universal or non-universal. 
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The total gaugino masses can be read as 

™ A™ \ — ™ tree _i_ ™ threshold , AMSB 
m a {mGXJT) — m \/2 + m l/2 + m a 

As will be seen, in G2 vacua the tree-level gaugino masses are suppressed relative to gravitino mass. 
This is due to the fact that the .F-terms of the moduli are suppressed compared to those of the matter 
fields and the gauge coupling functions are independent of the T a . This also happens in [6]. Because 
the tree-level contributions to m\/2 are smaller than 7713/2, one has to carefully consider the anomaly 
mediation and threshold corrections which could therefore be comparable to the tree level terms. In the 
following these contributions are evaluated at the unification scale. 

a. Tree-level contribution The universal (having assumed gauge unification) tree level contribution 
to the gaugino masses can be computed as [201 EI] 

</! = e K ^F%f vs , (87) 



where the tree-level visible sector gauge kinetic function is given in (82). As mentioned earlier, we assume 
a GUT group broken to the SM by a discrete choice of Wilson lines. Thus, at the unification scale the 
gauge kinetic function of 3 gauge groups in the SM are the same. This implies that the tree level gaugino 
masses are also unified. We find 



U - id h-n/P^ e V31 + 2-8(n/P)l-n/P 21 98 1 - n/P J ' 



h /I = — U.iOC —re v ' ' — — — — TTiq/o , 100 

l / 2 s l-n/P \31 + 2.8(n/P)l-n/P 21 98 1 -n/P J A/ y ' 

where we have dropped the gaugino phases (Nitj)/P. The result from analytical calculation is plotted in 
figure [5| As can be seen, there is a funnel at n/P ~ 9/14 where the F-terms vanishes at tree level. 

b. Threshold Correction The threshold corrections to the gauge kinetic function from the Kaluza- 
Klein modes are computed in [12] 

a^uT = Im(/ vs ) + 5 , (89) 



where 



*=f7L = f In (4 sin 2 (90) 
In the above, 71, is a topological invariant and q and u are integers determined by the topology of Xj. 



Then, we invert (89) to find 

Im(f vs ) = a^ lJT -5 = r l a^ lJT , (91) 

where now 

77 = 1 - «gut^ • (92) 

Summing up both contributions we fine 

m tree + threshold = 1 ^pig.^ (93) 

1 2zlm/vs 
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c. Anomaly Mediation Due to the substantial suppression of gaugino masses, it makes sense to 
take into account the anomaly mediated corrections [22]. These contributions are given by the following 
general formula at the unification scale 



rn 



anomaly 



«GUT 
47T 



(3C a - C a ) e K / 2 W + (C a - C«) e K l> £ F 1 Kj + 2 C a e^ 2 £ F I di In K 



vs 



(94) 



where C a = 0, 2, 3 and C a = ^ a = 33/5, 7, 6 are the quadratic Casimirs of the a th gauge group. In 
our case it can be written in a more compact form as 



m 



anomaly 



«GUT 
4-7T 



(3C a - C a ) + G 7 K 7 HS+c/l (c o - C a (l - 2 C/ 



where the parameter cj is defined as the following ratio 

<9/lnK vs 



ci 

At the minimum of the potential one obtains 

14 1 /5 15 1 

+ 



K 



HS+(7i 



(95) 



(96) 



8=1 

3 



0.46 



a=l 



-8/9n/P 
1-n/P 



+ 8.4(n/P) 



l + 2.8(n/P) 

4 7 
9 

4n/P 

9 ~ 91 + 2.4(n/P) 1-n/P 



n/P 
1 - n/P 



p e -AA(n/P) 



3 1-n/P V3 141 - n/P 

3 "/^ , 7 23.5(n/P) 2 _ /4 _ 7 1 3 - n/P 

1-n/P 9 1 + 2.8(n/P) ' Kl } \3 9 1 + 2.8(n/P) 1 - n/P 



(97) 



7 



1 



1 — n/P 



(98) 



Furthermore, one can see that for different moduli Sj and T a 



1 



l 



_ £a T a T a 

7 3V y 



1 + 



1 



2.8 7 (n/P) 



o | E a r n T a 3 + 8.4(n/P) 



7 /3 



7 /3 



i | - XgTaTa l + 2.8 7 (n/P) 



(99) 
(100) 



Putting all pieces together, one finds 

anomaly QGUT 



4tt 



b 3 x 



{3C a -C a )+Y,G i K l (C a -C a (l-2c i ))+Y,G a K a (C a -C a (l-2c a )) 



i=l 



a=l 



m 3/2 . (101) 



Here we have also dropped the phase factor which is identical to the one in the tree level contribution. 
The anomaly mediated contributions are numerically comparable to the tree-level terms but are not 
universal. Thus, the gaugino masses are non-universal. The anomaly contributions are also plotted in 
figure [5} 

In summary, the full gaugino mass at the unification scale depends on the parameters n/P, P, A, c, 
Im/ vs , 7 and r\ (or equivalently 5). The ratio of the gaugino masses to the gravitino mass just depends 
on Im/ vs , 7 and ij which are subject to the constraint (91). The total gaugino masses and the gravitino 
mass are plotted in figure [6] in a restricted range of parameters. 
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P=10, Ac~ 3/2 =100, a am -'=25, 17= 1, y=0.5 




P=10, AcT 3/i =100, q-gut - =25, ipl, r=0.5 



0.500 r 




0.001 



0.4 05 0.6 0.7 0.8 0.4 0.5 0.6 0.7 0.8 

n/P n/P 

FIG. 5. This plot shows gaugino to gravitino mass ratios. On the left we have separately plotted the tree level and 
AMSB contributions. On the right, they are summed up. It is clear that for all values of n/P where expansion 
makes sense, this ratio is smaller than one. 



2. Scalar Masses 



The general formula for the un-normalized scalar masses is [20 



m' 2 

m al3 



l m 3/2 + 



(102) 



I, J 



One has to canonically normalize the visible sector Kahler potential by introducing unitary normalization 
matrices. The scalar masses generically have a flavor diagonal contribution of order the gravitino mass 
and a flavor non-diagonal contribution proportional to gravitino mass. The latter is hard to calculate 
and it is argued above and in [6] they are small compared to the former. In the limit of a flavor diagonal 
Kahler metric and tuned cosmological costant, the masses for the canonically normalized scalars are 



nr 



a/3 



m 3/2 



Y,e KX+Ul F I F J d I dj\nK vs \6, 



a/3 ■ 



(103) 



In the minimum of the potential, one obtains the following expression in the leading order 



m 



I +^(l8.2(n/P) + 3 n,P 



6.5(n/P) 3 -n/P 18.2(n/P) 2 



I -n/P l + 2.8(n/P) I- n/P 1 + 2.8(n/P) 



1/2 



m 



3/2 



(104) 



The anomaly contributions to the scalar masses are suppressed relative to the order gravitino mass terms 
and we will thus neglect such contributions. In figure [7] the scalar masses are plotted against n/P and 7. 



3. Trilinear Couplings 



The un-normalized trilinear couplings are given by the following general expression [20 

w 



4/ _ p K/2 <Z V F 1 



(105) 
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FIG. 6. Masses versus n/P (upper left). Masses for fixed n/P versus 7 (upper right), n (lower left) and «gut 
(lower right). 
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FIG. 7. Scalar mass to gravitino mass ratio versus n/P (left) and 7 (right). 



In the case of a diagonal Kahler metric in the visible sector and under assumptions we made in scalar 
masses, the normalized trilinear couplings are proportional to the physical Yukawa couplings Y a ^ as 

A aM = Y a p y e KX+Ul / 2 WY, Gl [ Kl + dl lnY ^ ~ dl ln {K* S Kj S Ky S )] . (106) 
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As mentioned earlier, Y' a a are the un-normalized Yukawa couplings which arise from Euclidean mem- 
branes wrapping associative 3-cycles and connect isolated conical singularities. Due to the approximate 
flavor- diagonal nature of the Kahler potential, one finds 



Y 



_ K -(HS+ir 1 )/2 . , , ~ vs ~ vs ~ VSN-1/2 _ ^(hs+c/ i) /2 



(107) 

Having done some algebra, one can show that, at the minimum of the potential, the trilinear coupling 
at tree-level can be read as 



— Y a pj 



m 3 /2 



(108) 



The first two terms are the same as in (97) and the last term is the following 

n/P 



4 



a/37 



i=i 



l-n/P 



«HS • P 



1 



15/14 



14/3 

1 - n/P 1 + 2.8(n/P) ~ 3 ~ l-n/P 



V- 



73; 



(109) 



where ^73 is the volume of the 3-cycle connecting co-dimension 7 singularities where the chiral multiplets 
are localized. Generically, the trilinear couplings are the same order of gravitino mass. 

4- [i and B[i Terms 

Having set the bilinear term in the superpotential to zero, the general expressions for the normalized 
H and Bji are given by [20] 



(m 3/2 Z -e K ^Fi dl z)(K™ K V H: 



'VS z^VS\-l/ 2 
Id) ' 



(110) 



Bfi 



{2m 2 3/2 + V )Z- m 3/2 e K / 2 ^ F'djZ + m 3/2 e K / 2 £ F 1 (d z Z - Zdj In {K™K™ 



e K J- F T F J (djdjZ - djZdj In {KllKll 



VS &<-VS\-l/ 2 



1111 



The Higgs bilinear Z in the Kahler potential is treated in detail in {24j. In our study we will parametrize 
the \x and Bfi terms as follows 



fj, = Z^m 3 /2, 
Bfi = Z B ^m 



3/2 



(112) 
(113) 



where Z 11 and Z Bfl are 0(1) parameters. 



E. Physics at the Electroweak Scale 



In the previous subsection, all of the soft SUSY breaking parameters were computed at the unification 
scale. In order to connect to low energy physics and get the MSSM spectrum at the electroweak scale 
we need to RG evolve the masses and couplings from the unification scale down to the electroweak scale. 
The detailed discussion of renormalizing the masses and couplings down to electroweak scale is given in 
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Due to the anomaly mediated contribution, the gaugino masses are sensitive to the value of «gut- 
However, the value of ogut is only determined once we know the exact spectrum and run the gauge 
coupling up to the GUT scale. 

Although at the unification scale squarks and sleptons have approximately flavor universal masses, 
it turns out that RG evolution down to the electroweak scale increases the masses of the 1 st and 2 nd 
generation squark and sleptons and decreases the mass of the 3 rd generation. The former increase is 
mostly proportional to the gaugino masses and thus the effect is small. The latter decrease is due to 
large Yukawa couplings and is caused by the trilinear parameters. In particular, the right handed stop 
is typically the lightest of all sfermions. However, the [i parameter which determines the masses of the 
Higgsinos does not change much due to the standard non-renormalization theorems. 

As in [6] , it is convenient to work with an effective theory with sfermions and Higgsinos integrated out 
at the decoupling scale mj ec ~ 0(m 3/ / 2 )- The low energy effective theory below this scale is the Standard 
Model plus gauginos. To take into account the threshold effects of integrated out fields, we exploit the 
match and run procedure. In the rest of this subsection, we study the evolution of gaugino masses. 

The gaugino masses at the electroweak scale at one-loop level can be related to those at the unification 
scale by a RG evolution factor 

m a (m E w) = K a m a (m G UT) • (114) 

Whereas the running is done with tree-level matching at the decoupling scale, the RG evolution factor 
can be read as 



a a (m de 
A„ = 



«GUT 



( M^ M+g /^ M+g 
Qq(wiEw) 



(115) 



The constants ba M+g = 0, —6, —9 and 6a M+ ^ = 41/10, —11/6, —5 are the 1-loop /3 function coefficients of 
the gaugino masses and gauge coupling respectively. At the decoupling scale, the gauge coupling is given 
by 



a a (m dcc ) = a GUT + — In (— — J, (116) 
and at the electroweak scale it is 



a a 1 ( m -Ew) = a 1 ( m dec) + -^ hi( — ). (117) 

2ir \ rriEW ' 

are the beta function coefficients for the model above the decoupling scale e.g. for the MSSM with 
cxgut ~ 1/25, 6^ SSM = 66/10,1,-3. Combining equations (116) and (117) in order to remove gauge 
coupling in the decoupling scale a(md e c) we obtain 

-l, \ -l , ba M+g /m dec \ bf> M /m GUT \ . . 

a v ; GU1 2?r Vm E w' 2vr V m dec / v ' 

In figure [8] the 6^ SM are plotted for different values of m& cc and ogut- 

Thus, the RGE factors can be determined for a given unified gauge coupling acur and the decoupling 
scale mdec 



K a 



rSM+g /, SM+g 
i.SM+g / \-\o a B b a 

l-^a a (m EW )ln 



l + f^acUTln fe^ 



(119) 
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a OVI ~ 1 =20 (clashed), a avT ^=25 (solid), a am ~ l =30 (dotted) 
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FIG. 8. The /3 function coefficients for gauge couplings versus decoupling scale (left) and unified gauge coupling 
(right). 
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FIG. 9. The RGE factors for gauge couplings versus decoupling scale (left) and unified gauge coupling (right) for 
the MSSM. 



They are plotted in figure [9] for different values of m^ ec and ogut- 

Moreover, weak threshold corrections must be added in order to obtain the gaugino pole masses. 
Since Higgsinos are as massive as gravitino, there is a substantial threshold contribution from the Higgs- 
Higgsino loops which must be taken into account when computing bino and wino masses. Under some 
circumstances these corrections can be large enough to alter the nature of the LSP. The finite threshold 
correction is computed and given by 



Am finitc „ Q=l,2 



(120) 



This quantum correction can shift bino and wino masses up and down depending on the sign of fi. 
Furthermore, the electroweak threshold corrections from gaugino-gauge boson loops must be included 
especially for the case of gluino 



Atti-3 = —r- - ( 3 In ' 
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EW 



+ 5 m 3 . 



:i2i) 
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Taking all the various contributions into account, the gaugino masses at the electroweak scale are 



m a (m EW ) = ( r?m tree - m£ MSB ( 7 ))Ka + Am hnite (z e// ) . 



, finite i 



(122) 



In figures 10 and 11 gaugino masses at the electroweak scale are plotted. Evidently, the wino is the LSP 



in these vacua. 




IV. DISCUSSION 



We have considered a variation on the models considered in [B] in which a cubic superpotential for 
matter plays a significant role in moduli stabilization. This has allowed for a much wider range of gravitino 
masses than found in [6] , though for the phenomenologically viable window of 30TeV < m 3 / 2 < 500TeV 
the mass spectrum turns out to be practically identical to the models in [6]. In particular, we find that 
the neutral wino is the lightest supersymmetric particle. Note that a bino or Higgsino LSP would be 
inconsistent phenomenologically because their annihilation cross-sections are not large enough when they 
are produced non-thermally from the moduli decays before BBN [28]. The LHC signals of these models 
include multi top quark production from gluino pair production and decay as well as short track stubs 
from the chargino production and decay. 

The models under consideration show that supergravity interactions can generate vacuum expectation 
values for charged matter fields and thus break hidden sector gauge symmetries at a high scale, though 
the large vacuum expectation values which arise here really requires an understanding of higher order 
corrections in the matter fields which we cannot compute at present and leave for future study. It would 
also be interesting to consider different matter content in the observable sector e.g. Eq representations 
or the interesting single Higgs doublet models considered recently |29j . 
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